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ABSTRACT: 



This paper generalizes the classical momentum theory as usually applied 
to propellers, windmills and lifting rotors into a single unified treatment. 
It also extends the analysis to include the regime in which flow through 
part or all of the actuator disc is reversed with respect to the remote flow 
field. Dimensional analysis is used in a systematic manner to reduce the 
final results to their simplest and most significant forms. It is shown 
that the performance of these devices can always be represented by a single 
parameter family of curves in which the parameter expresses the extent to 
which the performance of the actual device falls short of the theoretical 
limit. Detailed algebraic solutions are derived in closed form for the 
performance of both the idealized zero torque and finite torque actuators; 
these represent performance limits which any comparable real device may 
approach but never exceed. A qualitative analysis is presented concerning 
the development of the vortex ring state. The fundamental dynamic stability 
of the ideal lifting rotor is also analyzed. 

While this paper deals with a classical topic in fundamental fluid 
mechanics, the approach employed is original and many of the results derived 
here are in a form that is essentially new. In particular, the solutions 
obtained for the reversed and transition flow regimes and for the idealized 
rotating actuator are believed to represent new formulations of these basic 
phenomena. 
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1. Objectives of the Analysis 

The simple actuator disc has long served as a useful idealization and 
standard of comparison for the performance, within a certain restricted 
operating range, of a broad class of fluid-mechanical devices. These include 
the ordinary aircraft propeller, the lifting rotor such as the helicopter 
rotor under certain conditions of vertical ascent or descent, and the 
unshrouded axial flow turbine or "windmill" as it is termed in this paper. 

The fundamental analysis presented here is applicable to all three of these 
devices which we designate here collectively as "actuators." 

The chief purpose of the present analysis is to extend the elementary and 
classical momentum theory of the idealized actuator so as the obtain theoretical 
performance limits over the entire possible operating range, including the 
regime in which the flow through part or all of the actuator disc is reversed 
with respect to the remote flow field I In all cases considered in this paper, 
however, the axis of the actuator remains parallel to the undisturbed flow far 
away. 

Another purpose of this paper is to utilize dimensional analysis in a 
systematic manner to reduce the final results to their simplest and most 
significant forms. It will be shown, in particular, that the performance of 
any of the foregoing devices can always be represented by a single parameter 
family of curves ; this parameter expresses the extent to which performance of 
the actual device falls short of the limit theoretically attainable from the 
comparable ideal actuator. 

This paper also considers certain fundamental aspects of the dynamic and 
hydrodynamic stability of the ideal actuator. 



1 



2 . Summary of Principal Concepts, Equations and Generalized Performance 
Charts 

The main results of this analysis are summarized in subsections 2.1 through 
2.11 inclusive. The detailed derivations are given elsewhere in this paper, 
but enough information is included in these subsections to define all results 
without the need to refer to the actual derivations. 

2.1 Nomenclature for Section 2 

All symbols which occur in the sximmary of principal results, subsections 

2.2 through 2.11 inclusive, are defined in the text when first mentioned and 
also, for ease of reference, in the summary below. Symbols are listed below 
according to the subsection and in the approximate order in which they are 
introduced in the text. A dash — after the definition signifies that the 
quantity is dimensionless. 

Subsection 2.2 

F = force exerted by actuator upon fluid, positive if in the same sense 
as the velocity of the undisturbed fluid. Ibf 
L = reaction exerted by the fluid upon a lifting rotor, always upward and 
positive by definition. Ibf 

V = axial velocity of undisturbed fluid relative to actuator disc, always 
non-negative by definition, ft/sec 

H = rate of vertical climb of a lifting rotor, positive upward, ft /sec 
T = sign changing operator defined by Eq. (2.2.3). -- 
Subsection 2.3 

O 

p = fluid density (constant), slugs/ft'^ 

D = diameter of actuator, ft 

= shaft power of an ideal actuator, positive if supplied to actuator, 
ft Ibf /sec 

P = actual shaft power, positive if supplied to actuator, ft Ibf /sec 
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= power loss due to all non-ideal effects, always non-negative. 



ft Ibf/sec 

N = shaft rotational speed, positive by definition, rev/sec 

O 

|jl = fluid viscosity. Ibf sec/ft 
Subsection 2.4 

Note: All dimensionless coefficients listed immediately below are defined 

mathematically in Table 2.4.1. 

Cp = conventional force coefficient in p, N, D system. 

Cp = conventional power coefficient in p, N, D system. 

J = conventional advance ratio in p, N, D system. 

Cp = power loss coefficient in p, N, D system. 
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Subsection 2.^ 

= axial velocity at actuator disc, positive if in the same sense as the 
undisturbed flow, ft /sec 

<P> T, o^, CT = sign changing operators defined by Eqs. (2.5). 

Subsection 2.6 

P, W, T, R = symbols denoting propulsive, windmill, transition and 
reversed flow operating modes, respectively. 

M = symbol denoting operating point which corresponds to generation of 
maximum possible windmill power lc^p|. 

B = symbol denoting operating point which corresponds to maximum drag 
force at zero net shaft power. 

a,m = parameters occ\irring in the cubic solution for thrust coefficient 



as defined in Eq. (2.6.4) of Table 2.6.2. 
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Subsection 2.7 



TR, IW = symbols denoting subregions of the transition mode corresponding, 
respectively, to positive and negative net shaft power. 

O', m = parameters occurring in the cubic solution for thrust coefficient 
as defined in Eq. (2.7.7) of Table 2.7.4. 

6 = parameter occurring in the cubic solution for velocity coefficient 
as defined in Eq. (2.7.7) of Table 2.7.4. 

Subsection 2.8 

5 = symbol denoting operating point which corresponds to static conditions 

(V = 0). 

6 = parameter occurring in the cubic solution for rate 

power as defined in Note 2 of Table 2.8.1. 

Subsection 2.9 

r^ = radius of non-swirling core in remote slipstream. 

cA 5 = maximum values of conventional thrust 

7 MAX 7MAX 

that can be attained at any specified 
an idealized finite torque actuator. 

R = slipstream radius, ft 
r = slipstream radial coordinate, ft 
0) = shaft angular velocity, rad/sec 

V = slipstream axial velocity, ft /sec 
o 

= slipstream peripheral velocity, ft/sec 
r = slipstream vortex strength outside core region, ft'^/sec 
H= slipstream advance ratio defined by Eq, (2.9-2) of Table 2. 9.1. 

A = slipstream loading parameter defined by Eq. (2.9-3) of Table 2.9-1. 

V = slipstream swirl parameter defined by Eq. (2.9-4) of Table 2.9-1- 
= relative radius of non-swirling core as defined by Eq. (2-9-9) of 

Table 2.9.1. — 

Q = auxiliary function defined by Eq. (2.9.7)' of Table 2.9.1. 



of descent at zero 
ft 

and power coefficients 
advance ratio J by 
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Subsection 2.10 



6L = perturbation in lift. Ibf 

W = total weight of rotorcraft. Ibf 

t = time coordinate, sec 

9 = dimensionless time coordinate defined by Eq. (2.10.4) of Table 

2 . 10 . 1 . -- 

C^p, = dimensionless performance coefficients in p, D, W 

system as defined by Eqs. (2.10.1) through (2.10.3) of 
Table 2.10.1. -- 

X, Y = dimensionless stability derivatives defined by Eqs (2.10.5) through 
(2.10.8) of Table 2.10.1. — 

p,, V = perturbation acceleration and velocity response constants, respec- 
tively, of control system as defined by Eq. (2.10.9) of Table 2.10.1. 

3 = exponential response rate of lifting rotor to a small perturbation. 
(Operation is stable if 3 is negative). 

2.2 The Lifting Rotor in Vertical Ascent or Descent 

In dealing with the lifting rotor in vertical ascent or descent, the 
rotor axis is placed vertically, and the sense of the flow through the actuator 
is so chosen that the reaction force F exerted by the fluid upon the actuator 
itself, whether it be a propulsive thrust, a windmilling drag or a power-on 
reversed thrust, always acts upward and thereby constitutes a useful lift L. 
These circumstances in no way affect the flow pattern itself, nor the basic 
method of analysis . 

It is convenient, however, to introduce two minor changes of variable 
when dealing specifically with the lifting rotor. In place of the usual axial 
force F which may be either positive or negative, we introduce the lift L which 
is always directed upward, and which we therefore take as positive by definition. 
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In place of the remote velocity V which is always non-negative by definition, 
we introduce the rate of climb H which is positive for vertical ascent, negative 
for vertical descent. (The case of slanting ascent or descent falls outside 
the scope of this analysis.) 

The required changes of variable are simple accomplished by the substitutions 

F = T L (2.2.1) 

V = T H (2.2.2) 

where t is a sign changing operator defined as follows, namely. 



T 




(2.2.3) 



By means of these substitutions, any result originally derived for the 
propeller or windmill can be converted to the corresponding result for the 
lifting rotor in vertical ascent or descent. 

2.3 Beal versus Ideal Actuators 

The present analysis is based on comparing the performance of any real 
propeller, windmill or lifting rotor with that of a corresponding idealized 
zero torque actuator. As illustrated in Section 11 of this paper, the above 
comparison can be made in more than one way. However, the scheme outlined 
below appears to be by far the simplest and most versatile. Hence this is 
the basic method which is consistently applied throughout this paper, except 
for the auxiliary discussion in section 11. 

In this method, the idealized actuator is arbitrarily assigned values of 
p, D, V and F which are respectively equal to those of the real device, by 
definition. This then suffices to fix the theoretical shaft power of 
the ideal actuator according to a definite relation which is derived elsewhere 
in this paper and which can here be symbolized in the form 
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(2.3.1) 



P. = P. (p, D, V, F) 



On this basis, the shaft power P of the real actuator, whether positive 
or negative, is always algebraically greater than the ideal shaft power P^. 



The difference between these two quantities, which we denote by symbol 
P^, represents a relative loss of power due to all non-ideal effects, including 
those associated with slipstream rotation, viscosity, nvunber of blades, detailed 



in general a complex function not only of p, D, V and F, but also of shaft 
speed N, of viscosity p,, and of various geometric parameters. It is necessarily 
non-negative. Thus symbolically. 



Therefore the shaft power P of the real actuator may finally be represented 
in the form 



Since the principal aim of this paper is to establish a simple performance 
limit which any real actuator can approach but never exceed, attention is- 
focussed primarily on the elementary analytical relation which is symbolized 
by Eq. (2.3.1) and which represents the idealized zero torque actuator itself. 

The more complex and obscure relation symbolized by Eq. (2.3.2), which relates 
the real actuator to the ideal, is defined conceptually in this paper, but is 
only partially analyzed. In a later section the overall power loss P^ is 
broken down ftirther into a rotation loss P^^ plus a residual loss P^. Of course, 
the latter quantity is in principle amenable to f\irther analysis, but this 
requires resort to methods which are either semi-empirical in character , or 



blade geometry, flow \insteadiness , separation effects and so on. Hence P^ is 



Pf = Pf (p, V, D, F, N, n, ) 2: 0 



( 2 . 3 . 2 ) 




(2.3.3) 
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else far more complex than those considered in this paper. 

In general, for any well designed actuator operating at or near its 
design point, the overall power loss should be comparatively small and well 
suited to estimation in any specific instance by accepted engineering methods 
based on the utilization of appropriate test data. 

Attention is also invited to the role of shaft speed N in the situation. 
Note that power loss is a function of N, as indicated in Eq. (2.3.2), but 
that the basic relation (2.3.I) for the ideal power P^ itself is actually 
independent of NI This fact is subsequently used to simplify the basic analysis 
in a very significant way. 

2.4 Dimensionless Performance Coefficients 

The analysis of any propeller, windmill or lifting rotor can be simplified 
and generalized by reducing all quantities and equations to non-dimensional form. 
Since the three fundamental dimensions of force, length and time are involved 
in this problem, it is necessary for this purpose to choose three of the basic 
parameters of the problem as dimensional reference parameters. This can be done 
in various ways. Four of the most important possibilities are summarized in 
Table 2.4.1. 

By far the most commonly used system of dimensionless coefficients is that 
labelled the "Conventional System" in the table. This uses p, N and D as 
reference parameters and results in the three familiar coefficients commonly 
known as the thrust coefficient C^,, power coefficient Cp and advance ratio J, 
as defined in the table. 

We observe, however, that while shaft speed N is a significant parameter 
of the overall problem, it is not involved at all in the fundamental relation 
of th^e idealized zero torque actuator as symbolized by Eq. (2.3 .I). Hence, 
although this relation can certainly be non-dimensionalized in. the conventional 
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Table 2.4.1 Summary of Dimensionless Performance Coefficients 
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p, N, D system, the results so obtained will not be expressed in terms of the 
fewest possible number of dimensionless parameters. In effect, the conventional 
system introduces an extraneous parameter N into the formulation of the basic 
relation (2. 3-1) and thereby creates one superfluous degree of freedom. 

The above complication is undesirable and can be avoided simply by dropping 
the conventional p, N, D system and adopting instead any one of the three alter- 
native schemes shown in Table 2.4.1. Notice that the various dimensionless 
coefficients that occur in these alternative systems can all be expressed very 
simply in terms of the familiar coefficients C_j C and J of the conventional 

r r 

system. These simple conversions are summarized in the table. 

The above three basic systems are theoretically equivalent in that they 
all express the same underlying relationship in terms of the same number of 
dimensionless parameters. The respective detailed results, however, are not 
necessarily equally apt and convenient in every possible application. As a 
rule, the specific context determines which one of these three systems expresses 
the pertinent results in the most appropriate form. 

It is also worth noting that under any circximstances in which parameter N 
reenters the problem in a significant way, the conventional coefficients C^, 

Cp and J once again become relevant. This is true, for example, of the idealized 
finite torque actuator as analyzed in detail in section 9> and summarized in 
subsection 2.9. 
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2.5 Sign Changing Operators 



Owing to the extremely wide range of operating conditions considered in 
this paper, the velocity through the actuator disc, the axial force F, the 
ideal shaft power and actual shaft power P are each subject to changes in 
algebraic sign. As an aid in dealing with these changes in an orderly fashion, 
we introduce a group of sign changing operators defined as follows: 



cp = 




+ 1 



(2.5.1) 



T 




(2.5.2) 



a. 

1 




(2.5.3) 



CT 




(2.5.U) 



The usefulness of these operators for generalizing the basic analytical 
relations will become apparent in the following subsections . 

2 .6 The Basic Solution in the p, D, V System 

The basic momentum and energy analysis shows that the complete operating 
envelope of an idealized zero torque actuator may be subdivided into four 
distinct ranges as indicated in Table 2.6.1. Each operating range has its own 
characteristic idealized flow pattern as shown schematically in Fig. 2.6.1. 

The two basic equations presented in Table 2.6.1 suffice to define the performance 
over this entire range. Provided only that the proper respective sign changing 
operators cp and be used as defined in the table, the single generalized 
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Table 2.6.1 Basic Actuator Performance Relations 



As Functions of Force Coefficient in p, D, V System 



Operating 

Mode 


Symbol 


Range 


Sign Changing 
Operators 


Basic 

Eq 




a. 

1 


Propulsive 


P 


0 ^ + CO 


+1 


+1 


(2.6.1) 


Windmill 


W 


- g (/3 - 1) ^ « 0 


+1 


-1 


(2.6.1) 


Transition 


T 


- ^ (»/~3 + l) ^ ^ ^ - l) 


Not 

Applicable 


(2.6.2) 


Reverse 


R 


- “ ^CiF ^ g + 1) 


-1 


+1 


(2.6.1) 



Basic Equations 



M odes P, W and R 



'IP 






Mode T 



n 



- C,„ + C, 



'IP 3/3 IF iPf 



( 2 . 6 . 2 ) 



Note: For definitions of coefficients, see Table 2.4.1. 
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ACTUATOR DISC 
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Fig. 2.6.1 SCHEMATIC DIAGRAMS OF OPERATING MODES. 




Eq. (2.5.1) applies to all three ranges designated as modes P, W and El The 
separate relation shown in Eq. (2.6.2) then applies over the remaining 
transition range, mode T. 

Eq. (2.6.1) is completely consistent with the usual elementary momentum 
analysis of the simple propeller or windmill. Hence it is not surprising that 
the same basic result is applicable to both. But, as Fig. 2.6.1 plainly shows, 
the range of fiilly reversed flow, mode E, involves a flow pattern which seems 
very much more complex than the simple patterns associated with modes P or W. 

It is therefore perhaps surprising that Eq. (2.6.1) can include all three of 
these modes within a single generalized expression. In addition to generalizing 
the classical theory in this striking way, the present analysis also yields the 
novel result for the transition range as represented by Eq. (2.6.2); this 
particular relation has no counterpart in the usual elementary theory I 

The basic overall relation defined jointly by Eqs. (2.6.1) and (2.6.2) is 
shown in Fig. 2.6.2 as a family of curves of C^p versus C^p, with C^p^ as para- 
meter. The curve for C^p^ = 0 represents the ideal actuator. This defines the 
limiting performance which any real actuator can approach but never exceed. 

The loss parameter C^pp expresses quantitatively how much the performance 
of any real device falls short of the theoretical ideal. It also provides a 
convenient basis for comparing the relative performances of various real actuators. 

Two critical points on the ideal curve are worthy of note. Point M shows 
that the magnitude of the power coefficient jC^p | generated by any windmill can 
never exceed |^ = 0.233. Point B shows that the magnitude of the drag coefficient 
|Cpp I which can be generated without the expenditure of net shaft power can 
never exceed =s O.605. For a real actuator of specified constant C^pp, these 

magnitudes decrease to (|^ - C^pp) and ( ^ - ^1= - C^pp), respectively. 
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Fig. 2.6.2 BASIC ACTUATOR PERFORMANCE CURVES IN THE p, D, V SYSTEM. 



It should be mentioned that if we set = 0 in Eq. (2.6.1), this 

analytical expression can generate not only the major part of the ideal per- 
formance envelope shown in Fig. 2.6.2, but also certain extraneous solutions 
which if they were shown here would fall within the interior region of the chart. 
These details are in face shown later in Fig. 6.2. It is demonstrated elsewhere 
in this paper that one of these branches , although it formally satisfies the 
equation, represents an unstable solution which cannot be physically realized 
in the present context. The other extraneous branches are theoretically stable 
but represent performance which is less than optimal. These extraneous solutions 
are not considered further in this stimmary of principal results. 

There is one singular operating condition, however, which does warrant brief 
comment at this point. That is the so-called ring vortex state, represented in 
Fig. 2.6.2 by point V. This condition is of particular interest in connection 
with the lifting rotor because it is associated with a catastropic loss of lift! 
The ts^ical flow pattern of the vortex ring state is shown schematically in Fig. 
2 . 6 . 3 . This operating state is considered further in the detailed analysis. 

The original solution summarized in Table 2.6.1 is based on 0^^, as the 
independent variable; this choice reduces the analytical expressions to their 
simplest form. In practice, however, it is frequently required to take C^p 
as the independent variable. Fortunately, the basic equations can in fact be 
solved for as a function of even though Eq. (2.6.1) happens to be a 

cubic equation in the variable C^p . When extraneous branches of the solution 
are discarded, the final results of this process can be summarized as shown in 
Table 2.6.2. Note that the single equation (2.6.1) is now replaced by the two 
expressions (2.6.3) and (2.6.4) which are applicable over different ranges. 
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Fig. 2.6.3 SCHEMATIC DIAGRAM OF LIFTING ROTOR 
IN VORTEX RING STATE. 
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Table 2.6.2 Basic Actuator Performance Relations 



as Functions of Power Coefficient in p, D, V System 



Mode 


Symbol 


Range | cp Eq. 

1 


Propulsive 


P 


i 

0 ^ (Cpp - c^pp) ^ i (2.6.3) 

' i 


Windmill 


W 


■ if ^ ^^IP ■ ^iPf^ ^ ° ’ (2.6.4) 

' I 


Transition 


T 

i 




N.A. 


(2.6.5) 


Reverse 


1 

1 

1 

: 




-1 


(2.6.3) 


Basic Equations 



Modes P and R 




(2.6.3) 



Mode W 

"if = f <"iP - °iPf> °°° [f " ” t ] 

/ 2T TT 

- — (C^p - Cppp) , 0 ^ a ^ ^ , and where m has the values -1 

and +1, respectively, over the following sub-ranges: 



(Table continues on next page) 
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Table 2.6.2 



Cont 'd 
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However, the analytical solutions summarized in Tables 2.6.1 and 2.6.2 are 
entirely equivalent, and both generate precisely the same family of ciirves, 
namely, that family previously shown in Fig. 2.6.2. 

It should be appreciated that the form of solution based on the p, D, V 

dimensional system as illustrated in Fig. 2.6.2 can be useful in any situation 

in which V is both known and non-vanishing. The windmill is perhaps the best 

example, because no windmill can generate useful shaft power, nor even a useful 

drag force, unless V is non-zero. Conversely, this particular reference system 

is virtually useless for dealing with any situation in which performance ■under 

static conditions, V = 0, is a significant factor because under these conditions 

all four of the dimensionless coefficients C,„ and become infinite 

IF IP IPf IN 

and comparatively meaningless. The important cases of the propeller under static 
thrust or the lifting rotor in stationary hovering flight are examples of situa- 
tions which cannot be usefully treated within the framework of the p, D, V 
dimensional system. To deal with these problems we must employ the other 
alternative systems. 
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2.7 The Basic Solution in the p, D, |p 1 System 



In the previous p, D, V dimensional system, the parameters p and D are 
innately positive, and V is restricted to non-negative values by definition. 

We note, however, that in the present dimensional system, shaft power P itself 
may be either positive or negative, depending on the particular mode of operation 
considered. It is essential that all dimensional reference quantities be 
restricted to positive values only, hence we take |p| rather than P as the 
appropriate reference quantity. 

We also note that the overall transition range, mode T, must now be sub- 
divided into two sub-intervals corresponding, respectively, to positive and 
negative values of net shaft power P. We denote these as modes TR and TW, 
respectively. 

Upon changing from the p, D, V system to the p, D, jpj system, the original 
pair of fvindamental equations as previously given in Table 2.6.1 are thereby 
converted into the corresponding pair of equations now shown in Table 2.7*1. 

Of course the change involved is purely a matter of algebraic form; the under- 
lying physical relationships are still precisely the same as before. 

We observe that Eq. (2.7*1) of the table is linear in but cubic in 

Cgp, whereas Eq. (2.7*2) is linear in 0^^ but cubic in Each equation is 

here solved for the linear term thereby expressing the final relations in their 

simplest algebraic forms, as shown. However, it is sometimes desirable to write 

both these relations for in the form of an explicit function of 

2F 2V 

this case it becomes useful to distinguish between the power input and power 
output modes. The solutions for these two cases are summarized separately in 
Tables 2.7*2 and 2.7*3* 

Conversely, it is sometimes convenient to be able to solve all relations 
consistently for as a function of This requires the inversion of Eq. 

(2.7*2) which is summarized in Table 2.7*4. 
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Table 2.7.1 Basic Actuator Performance Relations 



in p, D, |p] System 






Mode 


Symbol 


Range 


o 


CP 


Eq.’ 


Power Input Pfodes - 


Propulsive 


p 


+ 


2 ^^■‘^2Pf^ 


^ Cjp 2 0 


+1 


+1 


( 2 . 7 .: 


Reverse 


R 




2 


^ C^p ^ f (,^+l)(l-C2pp)^ 


+1 


-1 


( 2 . 7 .: 


Trans ition 
(Power In) 


TR 


+ 






+1 


N.A. 


(2.7.2 


Power Output Modes - 


Windmill 


W 


0 > Cpp ^ ^ i^-1) (l-«^2Pf^^ 


-1 


+1 


(2.7.3 


Transition 
(Power Out) 


W 


+ 




5 + » 


-1 


N.A. 


(2.7.2 


Basic Equations - 


Modes P, R, W 


^ (1 - oC^p^) 2a ^2F /o ^ ] 

2V oC^p TTcp (1 - oC2pp) 




Modes TR and IW 


2 O' (l - oCpp^) 

^2F ^ ^2Y c 

3/3 ^2V 





Note: For definitions of coefficients, see Table 2.U.I. 
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Table 2.7.2 Basic Actuator Performance Relations in Terms of 



Velocity Coefficient for Power Input Modes in p, D, |p| System 



Mode 



Symbol 



Range 



«P 



Eq. 



Propulsive 



0 ^ C2y ^ “ 



+1 



(2.7.3) 



Reverse 



0 ^ c^v ^ ^ 



§Z. (i _ c ) 

2tt ^2Pf^ 



1/3 



(2.7.3) 



Transition 
(Power In) 



TR 



|I (1 - C2j,p] ^'3 i ^ . 



W.A. 



(2.7.4) 



Basic Equations 



Modes P and R 



= r (1 - 



'2F I IT 



C2pf) 



2ll/3 



1 + 



1 + 2ncp ^2V 

27 (1 - 



1 , "27 



27 (1 - Cgpf) 



1/3 



1/3 



(2.7.3) 



Mode TR 



JI- „ 2 



(1 - Cgp^) 



'2F ■ 3v^ 2V 



"2V 



(2.7.4) 



Note; For definitions of coefficients, see Table 2.4.1, 
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Table 2.7-3 Basic Actuator Performance Relations in Terms of 



Velocity Coefficient for Power Output ivbdes in p, D, jp] System 



Mode 



Symbol 



Range 



Eq. 



Windmill 



W 






1/3 



^ ^ » 



(2.7.5) 



Transition 
(Power Out) 



W 






- ^ ‘^2Pf ^ ^2V ^ 



(2.7.6) 



Basic Equations 



Mode W 



^2F 



-/■ 



f " "2Pf) 



a ,2ttx 

3 " " <y) 



(2.7.5) 



21 ^2Pf TT 

where cos oi = J tt ^ j 0 £ a ^ , and where m has the values -1 

2tt n ‘~ 

2V 

and +1, respectively, over the following sub-ranges: 



m = -1 






1/3 

^ ^ H- « 


1 m = +1 
1 
1 

1 









Mode W 



^2F 3^ -2V 



V ^ 2 ^ ^2Pf^ 



(2.7.6) 



"2V 



Note: For definitions of coefficients, see Table 2.U,1. 
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Table 2.7.^ Basic Actuator Performance Relation in Terms of 



Force Coefficient for Transition Modes in p, D, 1 p| System 



Mode 


Symbol 


Range 


Of 


Eq. 


Trans ition 
(Power In) 


TR 


- 


^ (v^ + 1) (1 - 


2F 


+1 


(2.7.7) 


Transition 
(Power Out) 


TW 


- 


- 1) (1 + C^pff 


2F 


-1 


(2.7.7) 


Basic Equation - 


Modes TR and TW 


„ 3 /tt" “ °^2Pf^ 

cos p = - V 


where 



Note: For definition of coefficients, see Table 2.U.I. 
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The form of the solution is shown graphically in Fig 2.7 •! for the power 
input modes and in Fig 2.7-2 for the power output modes. In each case the 
solution consists of a family of curves of versus C with the power loss 
coefficient as parameter. Of course the curves for = 0 once again 

represent a limit that any real actuator may approach but never exceed. 

The upper part of Fig. 2.7-1 represents a propeller in normal operation. 

The lower part represents a propeller under conditions of power-on reversed 
thrust. Fig. 2.7-1 represents a windmill generating useful power and undergoing 
a drag force. We can think of this case as depicting a windmill of fixed size 
generating fixed shaft power under various wind velocities . Note the existence 
of a limiting velocity below which the specified power cannot be attained. This 

f py pi /o 

^ (l + Cppp) J . Also at any given 

value of above the critical, there are of comrse two solutions corresponding 
to low and high drag, respectively. 
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FORCE COEFFICIENT 




Fig. 2.7.2 BASIC ACTUATOR PERFORMANCE CURVES FOR POWER 
OUT PUT MODES IN />, D, |P| SYSTEM. 
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2.8 The Basic Solution for the Lifting Rotor In the p, D, L System 

Straightforward algebraic conversion of the original pair of fundamental 
equations of Table 2.6.1 into the p, D, L dimensional system produces the 
corresponding pair of equations for the lifting rotor as shown in Table 2.8.1. 
These equations show power coefficient C^p as an explicite function of rate of 
climb coefficient , with power loss coefficient as parameter. This 

arrangement not only yields the simplest algebraic result, but is also the 
most convenient for various other purposes as well. Of course, if necessary, 
both of these relations can be readily inverted to give C„- as an explicit 

oh 



function of (C^p - ^ 3 pf)j however, this again entails solving a cubic 
equation and the detailed results obtainable in this way are not included here. 

The fundamental relation defined in Table 2.8.1 is shown graphically in 
Fig. 2.8.1 as a family of curves. Once again, the curve for C^pp = 0 represents 
a performance limit that any real lifting rotor can approach but never exceed. 

Of particular interest are the i^o critical points marked S and B. Point S 
shows that the power coefficient C^p required to maintain steady hovering can 



never be less 



than V — = 0 . 79^ • Point B shows that the rate of d escen t at zero net 
shaft power as measured by 1 never be smaller than = 1.286. The 

changes of these values with increasing loss coefficient C^p^ is also plainly shown. 

It. is noteworthy that these results represent an elementary solution in 
closed form over the entire operating range! Hitherto, elementary solutions 
based on the ordinary momentum analysis have been available over those portions 
of the range which correspond, in the present terminology, to modes P and W 
only. Moreover, the present results embody the useful parameter C^p^ which 
facilitates the comparison of any real lifting rotor with the ideal, or with 
some other real rotor. 
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Table 2.8.1 Basic Lifting Rotor Performance Relations in Terms 
of Rate of Climb Coefficient in p, D, L System 



Mode 


Symbol 


Range 




T 


a. 

1 


Eq. ■ 


Propulsive 


P 


+ ® s >: 0 


+1 


+1 


+1 


(2.8.1: 


Reverse 


R 


° * «3H ^ 


-1 


-1 


+1 


(2.8.1 


Transition 


T 




N.A. 


N.A. 


N.A. 


(2.8.2 


Windmill 


W 




+1 


-1 


-1 


(2.8.1 


Basic Equations - 


Modes P, R, W 


P (2.8.1 




Mode T 

"3P = * "3S^ - V =3Pf 



Note 1: At = 0, Eq. (2.8.1) gives C^p = 



Note 2: At C^p = 0, Eq. (2.8.2) gives 



.-M 



/6 , 2h\ 

n <3 * 



where cos 6 



. = . 7m 

’j V _ • 



H 

'3Pf’ 2 



~ ^ 6 s rr . In the limit = 0, this reduces to 



3H • TT 

Note 3: For definitions of coefficients, see Table 2.4.1. 



30 




Fig. 2.8.1 BASIC ACTUATOR PERFORMANCE CURVES FOR LIFTING 
ROTOR IN p, D, L SYSTEM. 
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2 .9 The Idealized Rotating Finite Torque Actuator 



The ideal actuator mainly considered in this paper is essentially a zero 
torque device. In section 9> however, a detailed analysis is presented of an 
actuator which, although ideal in all other respects, operates at finite torque 
and at finite shaft speed. Such a device will always operate somewhat less 
effectively than the ideal zero torque actuator, although at high shaft speed 
and low torque, the difference will be small. In dealing with the finite torque 
actuator, it turns out to be advantageous to revert to the use of the conventional 
performance coefficients 0^, Cp and J based on p, N and D as dimensional reference 
parameters. Moreover, the relations among Cp, Cp and J can be expressed paramet- 
rically in terms of certain slipstream parameters as defined in Table 2. 9.1. 

The analysis of section 9 shows that the torque and rotation effects lead 
to the creation of peripheral velocities in the slipstream leaving the actuator. 

A solution is found for these peripheral velocities which maintains the axial 
velocity distribution vinchanged from what it would be for the ideal zero torque 
case. In other words the axial velocities at the disc and in the remote slip- 
stream remain constant and \iniform as before. The slipstream is found to 
consist of two distinct regions, namely, a small inner core of radius r^ in 
which the flow remains purely axial, and an outer annulus in which the flow 
is of the free vortex type. 

The chief res\ilts of this analysis are summarized in Table 2,9.1. Notice 
that this analysis of the rotating actuator is restricted to the simple modes 
P and W only; the extension of the analysis to include all modes appears to be 
feasible but has not been undertaken in this paper. 

In principle, the equations of Table 2.9.1 suffice to fix a family of 
curves of C„ versus J, with C_ as the parameter whose value changes from one 

r Jr 
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Table 2.9.I Conventional Performance Coefficients of Idealized 



Finite Torque Actuator in Terms of Slipstream Parameters 



Dimensional Slipstream Parameters 

R = Slipstream radius 

r = Slipstream radial co-ordinate 

(U = Shaft angular velocity (rad/sec) 

V = Axial velocity of undisturbed flow 

V = Slipstream axial velocity (constant) 

O 

= Slipstream peripheral velocity 
V^r = r Slipstream vortex strength (Constant over the region 
r^ r ^ R.) 

Dimensionless Slipstream Parameters 

^ = K Slipstream advance ratio 

(Vg - V) 

= A Slipstream loading parameter 

r 

— ^ = Y Slipstream swirl parameter 



0 



— = Relative radius of non- swirling core 



Range of Loading Parameter 




Positive values of A correspond to mode P, negative values to mode W. 
modes are not considered in this particular analysis. 



(Table continues on next page) 



(2.9.1) 

(2.9.2) 

(2.9.3) 
(2.9.^) 

(2.9.5) 

(2.9.6) 

Other 
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Table 2 . 9.1 



Cont ' d 



Aipciliary Functions 



c 



/i - y 1 



2 (h + ^) A 



(h + ^) A 



(2.9.7) 



1 + -^ (h + 2 ) A + |-(h + ^) 






V = C (h + ^) A 




_2_ 

CY 



Parametric Solutions for Conventional Performance Coefficients 




(2.9.8) 

(2.9.9) 



(2.9.10) 

(2.9.11) 



(2.9.12) 



Static Conditions and Maximum Loading Limits 

Under static conditions (V = O), the above relations simplify as indicated 
in section 9 * 13 ? with results as shown in Fig. 2 . 9.1. 

Under maximum loading conditions, they simplify as indicated in section 
9.1^, with results as shown in Fig. 2 . 9 . 2 . 

/ 

(End of Table 2.9.I) 

3 ^ 



curve to the next. The calculations required for this purpose are straight- 
forward but tedious. At the time of writing, they still remain to be carried 
out. When completed, the resulting family of curves will define idealized 
performance limits which provide a suitable standard of comparison for the 
performance of any real single rotation actuator. By application of the 
rationale explained more fully in section 9> these results may also be used 
to evaluate the idealized rotation power loss for any prescribed operating 
condition. 

Under static conditions (V = O), the general solution siunmarized in Table 
2 . 9.1 reduces to the special form given in section 9*13- The corresponding 
result is shown in Fig. 2.9.1 of this section in the form of a curve of the 
quanties Cp and Cp/Cp ^ plotted versus Cp. The intercept at Cp = 0 corresponds to 
the limiting performance at zero torque, and the drop-off of this relative 
thrust parameter with increasing power coefficient Cp shows the effect of 
rotation losses. 

Another special case of theoretical interest is that corresponding to the 
maximum loading theoretically attainable, as analyzed in section 9-l^- The 
essential results are summarized here in Fig. 2.9*2 which shows and 

Vmax 



plotted as functions of J. Note that this absolute performance limit 

7MAX 



is deduced piirely from fundamental dynamical principles and without any reference 
whatever to any empirical data or to such complex and highly variable details 
of the flow field as local lift coefficients, separation effects or the like'. 
These latter effects can and do reduce the performance of any actual device 
well below the limits shown in Fig. 2.9-2. The essential point, however, is 
that there is no feasible way to raise the actual performance above these 
limits . 
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THRUST PARAMETER THRUST COEFFICIENT 
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FINITE TORQUE ACTUATOR. 
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It is also noteworthy that Table 2.9.1 summarizes the performance of the 
ideal finite torque actuator in a comparatively simple algebraic form. In 
general these results yield information that goes well beyond what can be 
deduced from the usual elementary analysis which simply neglects rotation effects 
altogether . 
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2.10 Stability of the Ideal Lifting Rotor 



In this subsection, we summarize the dynamic response of an ideal 
lifting rotor (P^ = O) to small perturbations from equilibrium during vertical 



ascent or descent. At the equilibrium state, lift L is equal to the rotor 
craft weight W, but lift is subject to small perturbations 6L about the equilib- 
rium value . 

It is advantageous in this application to choose p, D and W as dimensional 
reference parameters. The corresponding dimensionless coefficients of power 
C__, climb rate C_£j, and lift C. are then as defined in Eqs. (2.10.1) through 

O-T o-u oil 

( 2 . 10 . 3 ) of Table 2.10.1. It also simplifies matters to express time in terms 
of the dimensionless coordinate 9 as defined in Eq. (2.10.h). 

The inherent stability characteristics of the lifting rotor are expressed 



rized in Eqs. (2.10. 5 ) through (2.10.8). 

The power input to the rotor or, in descent, the power output from the 
rotor may either be held constant, or may be modvilated by a control system 
in response to sensed values of the perturbation acceleration and/or pert\irbation 
velocity, according to Eq. (2. 10. 9 ). Constants p, and v express these two 
components of the power response . 

The final dynamic response of the rotor craft to a small perturbation 
is then exponential in character as shown in Eq. (2.10.10). A stable response 
corresponds to a negative value of the exponential factor 3 . 

The solution for 3 is given by Eq. (2.10.11). This is seen to depend 
on the functions X and Y which describe the inherent stability characteristics 
of the rotor itself, and upon the constants p and v which characterize the power 




which must be evaluated along the 




control system, if any. 
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Table 2.10.1 Stability of Idealized Lifting Rotor 
in Vertical Ascent or Descent 
(Pf = 0) 



Dimensionless Coefficients 



C3P = (2.10.1) . I (2.10.3) 



1 

D y p 



^ ^ (2.10.2) e = (2.10.4) 



1 

D J p 



i_ ft 

gD y p 



Actuator Derivatives at Equilibrium (C-^. = l) 

ob 

Modes P, W and R 



ac 



X = 



3P _ 



9C 



■| + 



"h 

T 



3H 



CT. 

1 



2 / TTCp 



y = 



^3P ''h 

5 Cl - 2 



2 / tt CO 






]ybde T 



X = ^ i T 0/ - 1 



Be 



3H 






H 



acsp 

Y - z= _ P . 

ob 



(2.10.5) 



(2.10.6) 



(2.10.7) 



( 2 . 10 . 8 ) 



(Table continues on next page) 
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Table 2.10.1 



Cont 'd 



Control System Coefficients 



6C 



3P 






1 



6C 




U = perturbation acceleration response constant (if any) ^ 0 
V = perturbation velocity response constant (if any) s 0 



(2.10.9) 



Dynamic Response to a Perturbation 



^3H = i 'Sh 






0 



where 



& = - 



X + V 
Y ^ 



Algebraic Signs of Actuator Derivatives 



X > 0 for all values of C 



3H 



Y > 0 if C-A > - ^ 
3-b 



Y < 0 if C- • < - ^ 



Stability 

Range of stable response for |J. = v =0 



0 < 0 for - < C„A < + ® 

*/rr 3H 



Stability Augmentation - 

1. Increasing y. increase stable range 

2. Increasing v - 

a) does not affect stable range 

b) increases response rate (whether stable or unstable) 



( 2 . 10 . 10 ) 



( 2 . 10 . 11 ) 



( 2 . 10 . 12 ) 

(2.10.13) 



(2.10.14) 



(End of Table) 
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Calcialations reveal that the rotor derivatives X and Y are as shown in 

Fig. 2.10.1. This leads to the conclusion that with p, = v = 0, the system 

is inherently stable over the range shown in Eq. (2.10.14). Inspection of 

Eq. (2.10.11) then shows that the range of stable operation can be extended 

at will by increasing the acceleration control constant p. Increasing the 

velocity control constant v does not affect the stability range itself, but 

augments the response rate, whether stable or unstable. 

The exponential response rate p which fixes the degree of stability or 

instability is plotted versus climb rate G • for several cases of interest in 

oh 

Fig. 2.10.2. 

It is to be expected that the stability characteristics of any real lifting 
rotor will generally be considerably less favorable than the stability of the 
ideal rotor as summarized here. 
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Fig. 2.10.1 STABILITY DERIVATIVES X AND Y FOR IDEAL LIFTING ROTOR. 



MODE W ' // MODE T MODE R MODE P 




Fig. 2.10.2 EXPONENTIAL RESPONSE RATE OF IDEAL LIFTING ROTOR TO A SMALL 
PERTURBATION. 




2.11 Auxiliary Figures of Lferlt 



The purpose of this brief subsection is simply to call attention to a 
number of ainciliary quantities which do not occur in the main theoretical 
development but which nevertheless provide useful figures of merit in certain 
applications. These include the concepts of efficiency, power effectiveness 
and force effectiveness as defined in section 11. 

3 . Basic Assumptions 

The purpose of this analysis is not to predict as accurately as possible 
the performance of any particular real actuator under specified conditions. 

It is rather to establish as simple as possible a theoretical performance limit 
which any real actuator might approach but which none could ever exceed. In 
general, the performance of any actual device involves physical mechanisms so 
complex that their effects cannot be adequately assessed by theoretical methods 
alone; it is necessary to refer also to the resvilts of experiment. Fortunately, 
the determination of a siiBple performance limit is a much easier task. For 
this purpose it is not only permissible but mandatory to neglect all complex 
but secondary effects , and to concentrate attention exclusively upon the 
irreducible essentials of the phenemena under study. 

For this reason we treat the fluid as incompressible and inviscid, we 
treat the flow field as possessing polar symmetry about the axis and as being 
steady, and we neglect possible hydrodynamic instability and turbulence effects. 
As a result of passing through the actuator, the fluid ultimately undergoes a 
change of axial velocity A V. We assume that this change develops asymptotically 
within the slipstream as the fluid moves away from the actuator and that, once 
established, the slipstream persists indefinitely. In particvilar, for any of 
the flow patterns involving flow reversal, or involving co-axial streams 
moving at different velocities and possibly even in opposite senses, we assume 
that these reversed or co-axial streams persist indefinitely without breaking down. 



For any real actuator operating in a real fluid, of course, the slip stream 
eventually disappears owing to the effects of hydrodynamic instabilities, 
turbulence and viscous dissipation. These effects are particular Important 
for cases involving flow reversal or co-axial streams, especially if the streams 
move in opposite senses. Nevertheless, even for real flows, it is usually 
possible to locate upstream and downstream cross-sections which are sufficiently 
close to the actuator that the slipstream is still well defined and yet far 
enough away from the actuator that the pressure in the stream is substantially 
equal to ambient pressure. For all practical purposes, conditions at such cross- 
sections are comparable to conditions "at infinity" in the idealized flow pattern. 
Hence the neglect of slipstream dissipation effects is not as drastic an ideal- 
ization as it might seem at first glance. Moreover, the usefulness of this 
approach has long been established in connection with the simple modes P and 
W. In this paper we are extending the same method to the other modes R and T. 
Admittedly, these latter modes, especially mode T, are more complex; nevertheless, 
for the special purpose of establishing a simple performance limit, the method 
appears to be fully justified. 

The idealized analysis also assumes that the velocities are purely axial at 
all reference cross-sections. Of course, any real propeller or other rotating 
actuator always exerts some torque upon the fluid and therefore sets up corre- 
sponding peripheral velocity components in the slipstream. Creation of the 
kinetic energy associated with these components absorbs shaft power. The 
resulting reduction of pressure over much of the slipstream also adversely 
affects the net thrust. Hence the peripheral velocity components degrade 
performance and are therefore eliminated from the idealized model, except of 
course in section 9 where this matter is considered in detail. 
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Similarly, a real propeller or non-optimal actuator may create non-imiform 
axial velocities in the far slipstream. Thrust is proportional to the mean axial 
velocity. Power required is proportional to the mean square of the axial velocity. 
Clearly, the ratio of useful thrust delivered to power required will be a max- 
imum only if the axial velocity is uniform across the slipstream. Consequently, 
the analysis of the idealized actuator assiunes uniform axial velocities across 
the far slipstream and across all reference stations. Pressures are also 
taken as uniform across these stations. 

In the transition mode T which involves two oppositely moving co-axial 
streams, the velocities and pressures are taken as uniform across each stream 
separately at each station; however, the values of velocity and pressure may 
in general be different for the two separate streams at a given cross-section. 

In the classical momentiun analysis of the simple modes P and W, it is 
assumed implicitly that the net axial force between the fluid which passes 
through the actuator and the ambient flow is zero or at least that it is negli- 
gible compared with the axial force exerted at the actuator disc itself. The 
corresponding assumption is here extended to modes R and T. In other words, 
the net axial force between any two co-axial streams is assiimed to be negligible 
in comparison with the axial force exerted on either stream at the actuator 
disc itself. The only exception to this occurs in section 9 which deals with 
the finite torque actuator; in this case it is necessary to admit a net axial 
force across the interface between the nonrotating core and the free vortex 
flow region. 

4. Additional Nomenclature for Main Text 

All symbols used in section 2 have previously been defined in subsection 2 . 1 . 
All other symbols occurring in this paper and not previously defined in subsection 
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2.1 are summarized below. Again the definitions are arranged according to 
the section and in the approximate order in which they are introduced in the 
text . 

Section 3 

AV = (Vg - V) = net overall increase in axial velocity occurring across 
the actuator, positive if in the same sense as the velocity V of the 
undisturbed fluid, ft/sec 

Note: For the important special case of static thrust defined by the condition 

V = 0, the velocities V^, Vg and AV, and the force F exerted by the actuator 
upon the fluid are all in the same sense and all have the same sign which is 
normally taken as positive. 

Section 3 

0, u, d, s = symbols or subscripts denoting the following locations, 

respectively: far upstream side of actuator disc, downstream 

side of actuator disc, and slipstream far downstream of actuator. 

V, V^, V^, Vg = respective axial velocities at the above locations, with 

V = V , . ft/ sec 
u d 

Pj P 5 Pjj P <3 = respective pressures at the above locations, with p = p. 

U CL S S 

Ibf/ft^ 

E = total energy per unit mass of fluid upstream of actuator, ft Ibf/slug 
£E = change in total energy per unit mass of fluid produced by actuator, 
ft Ibf/slug 

m = mass flow rate through actuator disc, slugs /sec 

P„ = power lost to slipstream by ideal zero torque actuator, ft Ibf/sec 

O 

Section 6 

= ideal power coefficient in p, D, V system. — 

X = ^ = velocity change ratio. 
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U = symbol designating that branch of the formal solution which represents 
physically unstable operation. 

S, 0, M, N, U, V, A, B, C, S’ = symbols designating significant points on 

ideal operating envelope as shown in Fig. 6.2. 



Section 7 

X = fraction of actuator disc area operating in the windmill mode d\aring 
operation in the transition range. 



S 

A, 




= slope of ideal operating curve 



C = subscripts referring to operation in windmill and reversed modes, 
respectively, at points A and C of envelope. 



Section 8 



= ultimate axial velocity of propagation of vorticity generated at 
actuator disc, ft/sec 



Section 9 
Subsection 9.2 

P^’ = shaft power of ideal finite torque actuator, ft Ibf/sec 

= idealized rotation power loss, always non-negative, ft Ibf/sec 
= residual power loss, always non-negative, ft Ibf/sec 
Subsection 

p = pressure at edge of slipstream at radius R (equal to press\ 3 re of 

p 

undisturbed flow). Ibf/ft 

Pg = pressiire within slipstream at radius r < R (lower than pressure of 

2 

undisturbed flow). Ibf/ft 
= peripheral velocity in slipstream, ft/sec 
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Subsection 9*^ 



2 

Pqq = pressure in central non-swirling core region of slipstream. Ibf/ft 

T) = g- = slipstream relative radius. 

Subsection 9 •10 

^!AX’ ^MAX’ ^MAX “ of the quantities A, Q, y under conditions of 

maximum loading attainable in an idealized rotating 
actuator. 

Subsection 9.1'y 

= relative angle of flow at radius r^ on downstream side of actuator 
as defined in Eqs . (9-15-1) or (9'15-^)- rad 
P = relative angle of flow at relative radius T] in remote slipstream 
as defined by Eq. (9.15.5). rad 



T)^ = mean radius ratio as defined by Eq. (9.15.6). 



Section 11 

T) = overall efficiency of real actuator as defined by Eq. (ll.l.). 

= efficiency of ideal actuator as defined by Eq. (11.1.2). 

e^^ = rotation power effectiveness as defined by Eq. (11.2.1). 

e^ = residual power effectiveness as defined by Eq. (11.2.2). 

ep = overall power effectiveness as defined by Eqs. (11.2.3) or (11. 5). 

e_ = force effectiveness as defined by Eq. (11.6.2). 
r 

e,, = braking effectiveness as defined by Eq. (11.6.3). 

Jd 

C.. j = maximum value of C.. for an ideal actuator at a specified value 
MAX 



of C 



IP' 



'IF 



MAX 



= minimum value of C for an ideal actuator at a specified value 

Ir 



of C 



IP' 
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5. Momentum and Energy Analysis of Idealized Zero Torque Actuator in Modes 



P, W and R 

For all three of the flow configurations P, W and R considered in this 
section, there are four major reference stations within each stream, designated 
respectively by subscripts 0, u, d and s. The velocities, pressures and total 
flow energies at these stations are designated as shown in Table 5.1. 



Table 5.1 Velocity, Pressure and Energy at Principal Stations 



Station 


Location 


Velocity 


static 

Pressure 


Total 
Energy 
Per Unit 
Mass 


0 


Far upstream 


V 


P 


E 


u 


At actuator, 
upstream side 


II 


Pu 


E 


d 


At actuator, 
downstream side 


^d 


Pd 


E + 


s 


Slipstream 


Vg = V + AV 


P 


E + AE 



It should be noted that the terms upstream and downstream in Table 5*1 
and elsewhere in this discussion always refer to the effective flow through 
the actuator, not to the ambient flow field. Thus for reversed flow in mode 
R, the upstream station u is on the aft side of the actuator disci Also note 
that by reason of continuity, the axial velocity must be equal at stations u 
and d and that the single symbol is therefore used for both. 

Bernoulli's equation can be applied firstly between stations 0 and u and 
secondly between stations s and d. Thus for any single mode we thereby obtain 



E 





(5.1.1) 
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( 5 . 1 . 2 ) 



(V + AV)^ _ % ^ 



2 



E + ZE = — + 
P 



P 



d 

2 



Upon subtracting ( 5.1.1) from (5.1.2) and rearranging, we obtain the difference 
in total energy and in static pressure across the actuator disc in the form - 



The present analysis is restricted to unshrouded configurations only. 
Thus the whole thrust F is produced by the above static pressure difference 
acting on the actuator disc area. However, the sense of the thrust depends 
on whether the sense of the flow through the disc is positive or negative. 
We therefore express the result in the generalized form 



where the operator cp assvires that the force always has the correct algebraic 
sign. 

A second independent expression for the thrust can also be found by 
applying the momentum theorem between stations 0 and s . For this purpose we 
need the absolute mass flow rate through the actuator which can be expressed 
in the form 



p® = (p^ 




AV 



(5.2) 




(5.3) 






(5.4) 



Consequently the required thrust may be finally written as 




(5.5) 



Comparison of Eqs. (5.3) and (5.5) now reveals that algebraically 




(5.6) 
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For the idealized actuator, the shaft power goes entirely into 
changing the total energy of the fluid which passes through it. Thus the 
algebraic shaft power equals the product of the algebraic energy change 
ZE as expressed in Eq. (5-2) with the absolute mass flow rate m as given by 
Eq. (5.^). Therefore, using Eq. (5*6) to eliminate V^, we finally obtain 
for the power 

?i = ^ P t) ( 5 . 7 ) 

At this point the key relations applicable to modes P, W and R can be 
summarized as follows: 



and 




( 5 . 8 . 1 ) 

( 5 . 8 . 2 ) 



(5.9.1) 



(5.9.2) 



Eqs. ( 5 . 9 ) define the fundamental relation among the principal quantities 
p. A, V, F and P^ in terms of the auxiliary parameter AV. While it is not 
mandatory to do so, it is in some respects convenient to eliminate this 
secondary parameter AV between the two relations (5.9). The resiilt may be 
written 

F^ = ^ p P^ (P^ - FV) (5.10.1) 
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We note that this equation is linear in V, quadratic in and cubic in 
F. Hence explicit solutions for V or P^ may be readily obtained, but an 
explicit solution for F, while entirely feasible, turns out to be more involved. 
We shall not undertake the cubic solution for F at this point. However, solving 
Eq. ( 5 . 10 . 1 ) for V gives 



V = 




_2_ 

TTCp 




(5.10.2) 



Likewise, solving Eq. ( 5 .IO.I) for P^ we obtain two roots. Nevertheless, by 
use of the operator a., the solution may be expressed unambiguously as follows. 




It is shown in a later section that the formal analytical solution defined 
by Eqs. (5.9)j or t>y Eqs. ( 5 .IO.I) or ( 5 . 10. 2), consists of four distinct'branches . 
Three of these correspond respectively to modes P, W and R. The fourth branch 
is extraneous; it represents a theoretical solution, mode U, which is actually 
unstable for any unshrouded configuration. One advantage of the particular 
form of solution defined by Eq. (5. 10. 3) is that by reason of the selective 
character of the operator this extraneous branch is eliminated. 

It is also instructive to analyze the basic momentum energy relations with 
respect to axes which, instead of being fixed with respect to the actuator 
itself as assumed heretofor, are fixed instead in the remote fluid. In this 
axis system, the actuator is moving with velocity V. The reaction force of 
amount F exerted by the fluid upon the actuator therefore delivers a thrust 
power of amount FV. Of course, this may be positive or negative thrust power, 
depending on the sign of F. 
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In addition, the velocities at stations 0 and s with respect to the new 



axes are now zero and AV, respectively. The corresponding total energies per 

" v2- 



unit mass are ^ and 
P 



P 



(AV) 



respectively. Hence the net change in total 

(av)^ 



energy per unit mass is the difference between these two quantities or 
This is simply the kinetic energy per unit mass of the slipstream with respect 
to the new axes. Since the slipstream flow is being "generated" at the rate 
of m units of mass per unit time, the creation of this kinetic energy in the 



slipstream consumes power at the rate P = m ' ^ . Of course the power P^, 

imparted to the slipstream represents a loss of useful power which is basic 
and unavoidable, even for an actuator which is, by definition, completely 
idealized. Moreover, since the thrust power FV, the power Pg lost to the 
slipstream, and the shaft power P^ must satisfy an overall energy balance, we 
may write at once that 



2 

- Fvj = m s 0 (5.11.1) 

By utilizing Eqs. (5*^) and (5*5)> this relation can readily be amplified to 
reveal that 





2 2 



S 0 



( 5 . 11 . 2 ) 



This relation in terms of the quantity can also be obtained independently 
by dividing Eq. (5*9.2) by (5.9*l) and rearranging. This agreement provides 
further confirmation of Eqs. (5.11.1) and (5.11.2). It also shows that the 
power analyses in the two different axis systems are mutually consistent. 

It is significant that the separate quantities P^ and FV can be either 
positive or negative, but that the difference (P^ - FV) must always remain 
non-negative. As explained above, this quantity represents an irreducible 
minimum loss of power which is unavoidable even with an ideal actuator. 
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Naturally, for any real actuator operating with corresponding values of p, D, V 
and F, the actual power loss must he even greater than this basic minimum. 

By rearranging Eqs. (5.10.1) and (5.10.3) we can also find two alternative 
and equivalent expressions for the basic power loss in terms that are frequently 
more convenient than the original form (5.11.2). They are, respectively. 





(5.11.3) 



and 




(5.11.4) 



Once again, the sign of the radical in Eq. (5.11.4) is specified by the 
operator in such a way as to eliminate the physically extraneous branch of 
the solution. 

Notice that the above results, Eqs. (5.8.1) through (5.11.4) have been 
established in a form which does not involve division by V in any term. Hence 
the solution for the iit5)ortant special case of static thrust may be found 
simply by substituting V = 0 into these equations. In this case we also find 
that since thrust power vanishes, 

= Pg (5.12) 



It may be remarked in passing that a fan, unlike a propeller, windmill 
or lifting rotor, is utilized for the very purpose of creating a moving slip- 
stream. Hence the power P„ imparted to the slipstream, which must be regarded 
as a dead loss for these other three devices, becomes the principal useful 
output for the fan I Apart from this single change in purpose and viewpoint, 
the analysis of the idealized fan is exactly the same as that given here for 
these other actuators . 
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The fundamental result of this section is the basic actuator relation, 
Eq. (5.10.3). Alternatively, we may choose either of the equivalent forms 
(5.10.1) or (5.10.2). The same basic relation is also defined by the pair 
of parametric equations, ( 5 . 9 «l) and (5.9.2), in which AV plays the role of 
auxiliary parameter. 

On the other hand for the exceptional case of a fan, the basic relation 

between output P_ and input P. is summarized parametrically by the pair of 
b 1 

equations , ( 5 . 11 . 3 ) and ( 5 > 10 • 3 ) j in which F now plays the role of auxiliary 
parameter . 

The foregoing fundamental results , in any or a 3 JL of the above forms , may 
be further simplified by non-dimens ionalization. This process reduces the 
number of parameters finally involved by three. The results previously sum- 
marized in section 2.0 were developed in this way. 

Eq. (5.10.1) which is cubic in F, is not solved explicitly for F in this 
section. However, the cubic solution has been worked out in detail and pre- 
sented in connection with several of the non-dimensional forms svimmarized in 
section 2.0. 

6. Stable and Unstable Branches of the Formal Solution 

It is useful at this point to non-dimensionalize the foregoing results in 
terms of the dimensional reference parameters p, D, V. For this purpose we 
define the following dimensionless coefficients, namely. 



'IF 






P. 

1 



p v3 



(6.1.1) 



( 6 . 1 . 2 ) 



X = 



jE 

V 



(6.1.3) 
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In terms of these coefficients, the detailed formal solution represented 
by Eqs. ( 5.8.1) through (5.11.U) may be summarized in dimensionless form as 
follows . 




( 6 . 2 . 1 ) 




( 6 . 2 . 2 ) 



cp = 




(6.2.3) 









'IPr 



(6.3.1) 



(6.3.2) 



*^1F^ ^IPi ("^IPi ■ ^1 f) 



(6.4.1) 



'IPi 



l^’i 
- 2 “ ^i 




MJ + 2_ 3 

Tl ncp ^IF 



(6.4.2) 



The last three of these equations as given below represent alternative 
expressions for the dimensionless power loss to the slipstream. 



C 



2 



\ 




2 ■ ncp 



(6.5.1) 



(6.5.2) 
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(6.5.3) 



^IF ^ 



/V \ /v^ 

The velocity ratios \-^j and 1 



= fi 



y. I as expressed by Eqs. (6.2) are plotted 
as functions of parameter X in Fig. 6.1. The corresponding relation between 
and C^p^, as defined parametrically by Eqs. (6.3) or explicitly by Eq. 
(6.4.1), is plotted in Fig. 6.2. 

Study of these charts discloses the existence of four distinct branches of 
the formal solution which are correlated with the previously defined modes as 
follows : 



Branch SO - Mode P (Propvilsion) 

Branch OMN - Mode W (Windmill) 

Branch NUV - Mode U (Unstable) 

Branch VOS' - Itode R (Reversed Flow) 

To see that branch NUV represents a physically mstable mode, we take note 
of the following facts. For this branch, parameter X lies between -1 and -2. 
Eqs. (6.2) and Fig. 6.1 then show that the flow through the actuator is in the 
same sense as the ambient flow while the final flow in the remote slipstream 
is reversed with respect to the ambient flow. Thus the flow reversal takes 
place after the fluid passes through the actuator. Moreover, Eqs. ( 6 . 3 ) and 
Fig. 6.2 show that C,- and C.^^. are both negative in this range. In other 
words the actuator sustains a drag force and extracts useful power from the 
stream which subsequently reverses and departs in a sense opposite to that 
of the ambient flow. 

While the above behavior wovild theoretically be possible within the confines 
of a properly designed channel with physical restraining walls, it cannot be 
sustained in the unrestrained environment of an unshrouded actuator. The 



59 



2.0-r 



I I- 

-3.5 1-3.0 




MODE I 



Fig. 6.1 







/ 



/ 
A 



/ 



/ 






/ 



- 1 . 0 " 



- I. 5" 



- 2 . 0 " 



MODE T 



-2.5" 



> 4 ^ 



MODE W 



44 



MODE P 



VELOCITY RATIOS IN SLIPSTREAM AND AT ACTUATOR 
DISC AS FUNCTIONS OF VELOCITY CHANGE PARA- 
METER X. 
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Fig. 6.2 IDEAL PERFORMANCE ENVELOPE IN p, D, V SYSTEM. 



mechanism of vorticity transport which causes this flow pattern to break down 
is discussed in some detail in section 8. It suffices here merely to point; out 
that a flow can be completely reversed with an unshro\ided actuator only by the 
expenditure of shaft power. Actually, the flow m\ist be reversed before it passes 
through the actuator, and must therefore enter it from the aft side. This 
requires the creation of a sufficiently intense suction on the aft side to turn 
the ambient flow inward and around toward the disc. Since branch NUV does not 
satisfy these various requirements, it must be discarded as physically irrelevant 
for any unshrouded configuration. 

The reader can verify that the analytical solution in the specific form of 
Eq. (6.4.2) automatically eliminates the above extraneous solution through the 
selective action of the operator o^. 

We note in passing that the detailed solution which is non-dimens ionalized 
in this section in terms of parameters p, D, V can also be non-dimensionalized 
in terms of other reference parameters, for example, p, D, jpj or, in the case 
of a lifting rotor, p, D, L. While the results so obtained differ in algebraic 
form, they express relationships that are theoretically equivalent to those 
summarized above. The form based on p, D, V happens to be most convenient for 
the present phase of the analysis. 

7 . A New Solution for Transition Range, Mode T 

Since the transition from mode W to mode R cannot actually take place 
along the unstable branch NUV of the original formal solution, we are compelled 
to formulate an alternative hypothesis to apply over this range. We now propose 

the following novel yet simple solution for this case . 

TT 2 

Imagine the total disc area D to be subdivided into two concentric regions, 
namely, an inner circle of area x ^ and an outer annulus of area (l - x) ^ 
where the transition parameter x lies somewhere in the range between 0 and 1. 
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We assume that the flow over the inner region corresponds to some as yet 
unspecified operating point A of mode W, and that flow over the outer region 
corresponds to some as yet unspecified operating mode point C of mode R. 

Consider the following relations. 

For the inner flow - 



-1 ^ ^ 0 




For the outer flow - 



(7.1.1) 



(7.1.2) 



(7.1.3) 



- » ^ Xc ^ - 2 



( 7 . 2 . 1 ) 



^IFC 



TT 

IT 




( 7 . 2 . 2 ) 



^IPiC “ ■ t 




(7.2.3) 



The corresponding net overall values for the combined flow are found simply 
by adding the respective contributions from the inner and outer streams . The 
results are 



^IF “ ^ ^IFA + (1 - ^iFC 



(7.3.1) 



^IPi 



= X C 



IPiA 



+ (1 - x) C 



IPiC 



(7.3.2) 
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Upon eliminating x between these last two relations we obtain 



“^IPi ‘^IPiC 



‘^IPiA ~ ^IPiC 

p - P 
IFA ^IFC 




(7.U) 



The resiilt shows that for fixed values of and the relation between 

Ct and C,„ is linear over the transition range. 

IPi IF 

The problem that remains is to find the values of the parameters and 

X such as to optimize actuator performance over the transition range. This 
u 

means that for any given value of C,-,, the power coefficient whether positive 

-Lj ±r 

or negative, should be algebraically as small as possible. Inspection of Fig. 

6.2 shows that this requirement is satisfied simply by choosing as the transition 
line the straight line segment AC which is tangent to the original analytical 
curve at the two points A and C as shown. 

In this connection we note that the slope of the original curve at these 
two points can be expressed in the form 



^d C 



IPi' 



d C 



IF, 



= ^A = 






3K 



il + 



TTT^ 



(7.5.1) 




(7.5.2) 



Furthermore, the slope of the straight line segment joining points A and 
C can be reduced to the form 




(7.5.3) 
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Now the optimizing condition requires simply that 



s = = Sc 



(7.5.4) 



Simultaneous solution of Eqs . (7*5) gives the following important results - 



S = = -1 



( 7 . 6 . 1 ) 



\ = -2 ) = - 0.845 



( 7 . 6 . 2 ) 






(7.6.3) 



Therefore Eqs . (7.1) and (7.2) give at point A - 



°1TA = - f 



(7.7.1) 






(7.7.2) 



At point C 



- 5 - ^•‘*3 



(7.7.3) 



'^IPIC - " i 



(7.7.4) 



Consequently, the resulting solution over the transition range as represented 
by Eq. (7.4) can now he reduced to the final form 



'IPi 



TT 

3 ^ 



'IF 



( 7 . 8 . 1 ) 
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which applies over the range 



- ^ + 1) ^ g - 1) 



( 7 . 8 . 2 ) 



8. Propagation of Vorticity and Development of the Vortex Ring State 

Consider conditions at station s far downstream of the actuator. The slip- 
stream at this location is moving with uniform velocity V_ , the ambient flow 
with uniform axial velocity V. In the idealized inviscid case, the two irro- 
tational flow regions are separated by a vortex sheet of infinitesimal thick- 
ness in the form of a cylindrical surface. The vortex lines are circumferential 
the surface consists of a continuous distribution of vortex rings. According 
to the Helmholtz laws of vortex motion for an inviscid fluid, the circulation 
moves with the fluid element which exhibits it. These rotating fluid elements 
which constitute the interfacial vortex sheet themselves move with a velocity 

V which is the algebraic mean of the velocities V and V_ on either side of 
CU b 

the interface, thus 



In other words the velocity of propagation of the remote vorticity is the 
same as the flow velocity through the disc itselfl This applies also to 
the transition mode T if we consider the inner and outer streams individually. 

We may view the vortex rings as vorticity shed originally from the actuator 
disc and propagated at an eventual axial velocity 7”^^ = Normally, the pro- 

pagation is in the sense leading away from the actuator, which tends to produce 
a stable flow pattern. If the sense of propagation should become reversed for 



= |(V + Vg) 



( 8 . 1 . 1 ) 



By comparing this with Eq. (5-8.2) we note that 




( 8 . 1 . 2 ) 
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any reason so that the vortex rings were to move toward the actuator disc, 



the tendency would be to disrupt the original flow pattern. All of the 
vorticity wovild tend to collect and become fixed at the actuator, thus 
leading to the establishment of the vortex ring regime as shown earlier in 
the schematic diagram Fig. 2 . 6 . 3 * 

This idea can easily be reduced to explicit mathematical terms. Consider 
firstly flows like modes P and W in which the slipstream continues in the same 
sense as the undisturbed flow. In this case 



^=(1 

V '' 



+ x) > 0 



and stability in the above sense requires that 



( 8 . 2 . 1 ) 



= (l + >0 (8.2.2) 

Next consider flows like mode R in which the slipstream is reversed with respect 
to the undisturbed flow. In this case 






< 0 



and stability requires that 



( 8 . 2 . 3 ) 




( 8 . 2 . 4 ) 



The above conditions may be consolidated into the single statement that, 
for stability, it is necessary that 
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(8.3.1) 



It can be seen that the stability constraint (8.3.I) is satisfied for »n 
values of X except those lying within the range 

-2 ^ X ^ -1 (8.3.2) 



But this is precisely the range of the mode U which was previously identified 
as the unstable branch of the formal solution I 

Of course, the two end points identified in Eq. (8.3.2) are points of 
neutral stability. We identify these as points V and N, respectively. Note 
that 

At Point N X = -1 (8.U.1) 



C 



IF 




(8.4.2) 




(8.4.3) 



At Point V 



X = -2 



(8.5.1) 



C 



IF 



= 0 



(8.5.2) 



C 



IPi 



= 0 



(8.5.3) 



It is of interest to compare this latter point with another significant 
operating condition, namely, point 0. Thus 
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At Point 0 



X = 0 



( 8 . 6 . 1 ) 




( 8 . 6 . 2 ) 



C 



IPi 



= 0 



( 8 . 6 . 3 ) 



This brings out the very interesting fact that although points 0 and V 
are both located at the origin in the versus plane, they represent 

quite distinct operating conditions. In particiilar, point 0 is a state of 
positive stability while point V is a state of neutral stability. In fact, 
point V represents the vortex ring state previously mentioned. A rotorcraft 
which is inadvertantly brought to this operating point suddenly suffers a 
catas tropic loss of lift! 

Consider in this light the operation of an actuator, initially in mode 
R, at gradually decreasing levels of power coefficient Assume that when 

the actuator reaches point C , operation is not shifted along the optimim 
transition path CBA, but is controlled so as to continue along in mode R 
toward point V. It can be seen that as the operating point moves closer to 
point V, the margin of stability in the above sense becomes less until at 
point V itself, the point of incipient instability is finally reached. It is 
evident that for operation unduly close to this poiai, a quite small pertur- 
bation may suffice to disrupt the precarious stability, and lead to the 
vortex ring state. 

Moreover, if the attempt were nevertheless made to continue the above 
gradual changes further into the unstable region VUN, the resulting unstable 
flow pattern would inevitably collapse, to be succeeded in most cases by the 
vortex ring pattern. 
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In contrast with this, during a shift along the optimum transition path 
CBA, operation too near the points of incipient instability, V and N, is 
avoided. Hence path CBA is theoretically stable in this sense as well as 
being the transition path of highest performance. 

Of course, for a real actuator, there may be various other types of 
instabilities to consider. Fortunately, these are usually of a less drastic 
kind. Even for the idealized actuator, however, there are other types of 
instabilities as well. One particular type is analyzed in section 10 of 
this paper. In general, however, the phenomena of hydrodynamic instability 
ir real systems are extremely complex and cannot be analyzed by the elementary 
methods of this paper. 

Rotation Effects on Actuator Performance 
9-1 The Idealized Finite Torque Actuator 

The ideal actuator considered so far involves purely axial flow in the 
remote slipstream and is, in effect, a zero torque device. In principle, this 
condition could be approximated in various ways . One method would be to use a 
pair of co-axial counter-rotating rotors which exert equal and opposite torques 
upon the flow. Another method would be to use a rotor and a co-axial stator 
which exert equal and opposite torques. A third method would be to utilize 
a single rotor at rotational speeds so high that the necessary power could be 
transmitted at torque levels low enough to be negligible. 

In most cases, the above methods are for one reason or another impractical, 
and it becomes useful to consider in somewhat greater detail the optimum per- 
formance which can be attained by a single rotation actuator which transmits 
appreciable shaft power at limited rotational speed and therefore with signi- 
ficant torque effects . In general such a device , which we term here the 
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finite torque actuator, even if idealized in every other respect, will always 
perform somewhatless effectively than the comparable idealized zero torque 
actuator . 

9.2 Idealized Rotation Power Loss 

Consider three actuators all of which operate with the same values of p, 

D, V and F, respectively, but with different values of shaft power P. The 
first is the ideal zero torque actuator with ideal shaft power as previously 
analyzed. The second is an idealized finite torque actuator and the third is 
the real actuator; both of these latter' devices operate at the same shaft speed 
N. The shaft power Pj^ of the idealized finite torque actuator equals (P^ + P^^) 
where Pj^, which is always non-negative, is here termed the idealized rotation 
loss. The shaft power P of the real actuator was previously written as equal 
to + P^) where P^, which is always non-negative, denotes the total power 
loss. This same total shaft power P can now be represented in somewhat more 
detail by the expression (P. + P„ + P ) where P , which we term here the 
residual loss, is again a non-negative quantity. This analysis resolves the 
total power loss P^ into an idealized rotation loss P^ plus a residual loss 

P . 

r 

9.3 Conventional Coefficients in Terms of Slipstream Parameters 

In the preceding discussion of power losses, it is necessary to distinguish 
between various power components and in this context the total shaft power of 
the idealized finite torque actuator is denoted by the somewhat elaborate 
symbol P^'* However, in this and subsequent subsections, the discussion is 
restricted solely to the analysis of the idealized finite torque actuator itself, 
hence we elect to denote this same power P^' hereafter by the simpler and more 
convenient symbol P. Consequently, the significant parameters of our problem 
are p, D, V, F, N and P. Note that this list now includes N, which is 
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non- vanishing 5 finite and positive. Hence it becomes advantageous to revert 
to the use of the ‘conventional performance coefficients Cp and J based on 
p, N and D as dimensional reference parameters. 

It is shown in a later subsection that the performance of the idealized 
finite torque actuator can be analyzed in terms of conditions in the remote 
slipstream, station s. It is necessary, of course, to distinguish between 
the diameter of the slipstream, which we here denote as 2R, and the diameter 
D of the actuator disc itself. It is also convenient to express shaft speed 
both- in conventional units of rev/sec as denoted by N and in scientific units 
of rad/sec as denoted by o). With this notation, the required performance 
coefficients may be expressed as follows: 



= 



F 






' F ^ 

2 ^ 
,p u) R y 



(9.3.1) 



= 



p D' 



r(f) 



p 0)^ R^' 



(9.3.2) 






(9.3.3) 



9.4 Basic Slipstream Analysis 

An actuator which transmits shaft power P to an initially nonrotating 
stream, and which itself rotates at shaft angular velocity u), necessarily 
exerts a torque of amount P/tu upon the flow. Assigning polar symmetry for 
the idealized flow model, we conclude that the above torque sets up in the 
remote slipstream a distribution of tangential velocity which varies in 
some manner as a function of radius r. In general, the axial velocity Vg 
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in the remote slipstream can also vary as some function of r. The stream- 
surfaces of this flow are straight co-axial cylinders . 

The pressure at the outer edge of the slipstream, radius R, equals the 
ambient pressure p. For a rotating slipstream, the pressure p_ at any interior 
radius r is below ambient. Considerations of radial balance of forces show 
that the amount of pressure drop is 




(9.4.1) 



Applying the energy equation between a point in the slipstream and a 
point far upstream of the actuator, we find the net energy input to the flow 
to be 







(9.4.2) 



The energy input can also be written in terms of shaft power input per 
unit mass flow, where shaft power is expressed in terms of torque and angular 
velocity, and where torque is written in terms of change in angular momentum. 
The result is simply 



AE = u) r ( 9 . 4 . 3 ) 

Upon combining these last three equations, we obtain the basic relation 

R 

V ^ r V ^ 

= (u V r ^ + I dr (9.4.4) 

u 2 J T 

r 

This shows how the distributions of V and in the slipstream must be related 

U b 

for any assigned values of the parameters V and co. 
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9.5 Idealized Velocity Distribution in the Slipstream 



Our basic aim is to find a distribution of V which will leave Y„ 

u S 

unchanged from the result obtained for the original zero torque actuator. 

In other words, we require a distribution which produces a uniform axial 
velocity Vg across the slipstream. Moreover, it must be possible to prescribe 
the magnitude of Vg arbitrarily. 

It turns out that the required solution consists of two distinct regions. 
There is firstly a non-swirling inner core of radius r^. Outside the core, 
there is an annular region of simple free vortex flow. Thus 

For O^r^r^ V r = 0 (9.5-1) 

0 u 



For 



r_^r^R V r = r = constant (9.5-2) 

0 u 



abstituting from Eq. (9-5-2) into (9-^-^) gives the result 




(9-5-3) 



With all other quantities arbitrarily prescribed, this equation fixes the 
characteristic swirl constant F of the outer free vortex flow. 

The total pressure drop through the outer annulus can now be fovind by 
substituting the free vortex relation (9-5-2) into the radial balance equation 
(9-^.1)- The result is 




(9-5-4) 



For the flow through the nonrotating inner core, there is no net energy input 
and the energy equation gives 
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(9.5.5) 




By eliminating the pressure term between these last two expressions, we obtain 
a solution for the core radius in the form 



R 



= 1 + 



'0 




(9.5.6) 



where F is now known from the prior solution of Eq. (9.5.3). 

Eqs. ( 9 . 5 . 1 ) through (9.5.6) completely define all features of the slip- 
stream flow which correspond to arbitrarily specified values of the parameters 
Vg, V and lo. 

9.6 Slipstream Thrust Integral 

Applying the axial momentum law to the propeller, we can now write the 
net thrust in the form of an integral over the slipstream as follows; 



R 



r 



F = 2rrp 

W 

0 






r 

/ 



r dr 



(9.6.1) 



The first term inside the integrand represents the basic thrust produced 
by the momentum increase across the actuator. The second term, which can be 
recognized as the pressure drop integral of Eq. (9.^.1)? represents a loss 
of thrust associated with the lowered pressures near the axis. This equation 
also indicates why it is necessary to maintain a nonrotating core of adequate 
radius r^. For if r^ be allowed to decrease unduly, the unfavorable negative 
pressure integral eventually dominates the solution. The non- swirling core 
of radius r^ eliminates this difficulty. 
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In fact, over the core of the present solution, the integrand of the 



thrust equation (9-6. l) reduces as follows 




2 




2 



> 0 



( 9 . 6 . 2 ) 



It is interesting to observe from this resixlt that the non-swirling core still 
makes a net positive contribution to the thrust despite the fact that it 
receives no energy input and undergoes no pressure rise in passing through the 
actuator I The reason, of course, is that it discharges as an axial jet into 
the region of reduced pressure created by the outer swirling flow. It exerts 
this thrust not as a pressure i ease across the actuator, but as a suction 
force distributed along the convf.rging stream surface which separates the 
inner and outer flows . The outer flow then transmits this force to the 
actuator disc itself. 

9.7 Slipstream Power Integral 

By integrating the energy relation (9.^ *3) with respect to mass flow rate 
across the entire slipstream, the following expression is obtained for shaft 
power, namely, 

R 

P = 2tt p U)J Vg r^ dr (9.7.1) 

r 



Likewise, 



from ( 9 . 4 . 1 ) and (9.4.2) we can obtain the alternative form 




(9.7.2) 
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The form (9.7.I) is much simpler to compute , but (9.7.2) is interesting 
theoreticeO-ly, especially in conjunction with the force equation (9.6.I). 



Notice that the effect of the terms in V is to decrease the thrust delivered 

u 

and simultaneously to increase the shaft power required I The adverse effect 
of slipstream rotation on overall performance is at once evident from these 
formulas. Fortunately, at sufficiently high shaft speeds, these effects 
become small. 

9.8 Non-dimensional Slipstream Parameters 

In applying the foregoing results it is useful to non-dimensionalize all 
slipstream quantities and relations initially in terms of the three dimensional 
reference parameters p, u) and R. In this connection we introduce the following 
dimensionless quantities, namely, 



V 



Slipstream advgince 
ratio 



(9.8.1) 




Slipstream loading 
parameter 



(9.8.2) 



r 




Slipstream swirl 
parameter 



(9.8.3) 



uiR 




Slipstream 

relative 

radius 



(9.8.4) 



Integration of Eqs. (9.6.I) and (9. 7.1) in this system then yields 
solutions for the following two coefficients, namely. 



^ — jj- = Slipstream force 

p U) R coefficient 



(9.8.5) 



1 — - = Slipstream power 

p R^ coefficient 



(9.8.6) 
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9.9 Slipstream Contraction or Expansion 

A glance at Eqs . (9*3) now shows that to convert the slipstream force and 

power coefficients into the conventional force and power coefficients C_ 

and Cp, and to convert slipstream advance ratio h to the conventional advance 

/2r\ ^ 

ratio J, it is necessary to know the slipstream area contraction ratio I— I 
Applying the continuity condition between actuator disc and slipstream 
gives the general result 




d/2 



I 



V r dr 
d d d 




r dr 



(9.9.1) 



However, for the specific solution developed in this section, the addition 
of the peripheral velocity components in the slipstream does not affect the 
original axial velocity . In fact the solution was specifically tailored to 

O 

meet this requirement. For incompressible flows having polar symmetry, a free 
vortex motion superimposed on an arbitrary eixial motion does not affect the 
latter. This principle is well known in fluid mechanics and is simply utilized 
here without detailed derivation or proof. We merely note that the superimposed 
free vortex motion does not affect the axial velocity distribution V^. Also, 
since and V, both happen to be constants, Eq. (9.9.l)simplifies to 

D Q 





Vg V + AV 



(9.9.2) 



Upon non-dimens ionalizing this in the p, (u, r system we obtain finally 



, HI 

iDy H + A 



(9.9.3) 
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9-10 Swirl Intensity, Loading Limit and Net Operating Range 

Fiirther details of the non-dimensionalized solution now work out as 
follows. The basic relation (9* 5*3) which fixes slipstream swirl intensity 
takes the form 



H+^)A=Y-iY^ 



(9.10.1) 



It is useful to solve this for \ the quadratic formula, and to arrange 
the resulting solution in the following special form, namely. 



Y = C + ^] A 



( 9 . 10 . 2 ) 



where the auxiliary parameter Q is defined by the expression 



l-Jl-2 

& \ / A\ 



A 



H + I A 



(9.10.3) 



/ 



It is also convenient for some purposes to expand this solution into 
series form as follows. 



C = 1 + i (k + A + ^ |h + A^ + 



(9.10.4) 



Inspection of Eq. (9. 10. 3) discloses an important physical operating limit 
Since Q is always real, the quantity under the radical can never become negative 
We conclude that 



-I 



_ 1. 



whereupon 



and 



MAX 



^^MAX ^ 



v)max = ^ 



(9.10.5) 



(9.10.6) 



(9.10.7) 
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Solving Eq. (9. 10. 5) for the loading limit by the quadratic formula 

gives the result 



'max 



= + 



y 



1 + H - K 



(9.10.8) 



There is also a negative root for A which would be significant if we 
were considering reversed flows, mode R, in this analysis. However, the 
present discussion is limited to the simple modes P and W only, and in this 
context the negative root is not relevant. 

On the other hand there is another restriction on A arising from the very 
fact that the analysis of this section ^ specifically limited to modes P and 
W. In terms of the earlier parameter we may say that the present study is 
restricted to the range 

_2 ~ ^ ^ X ^ » (9.10.9) 

\ ^ / 

where 

X = = -^ (9.10.10) 



Upon combining these conditions , we find that the present analysis is 
restricted to the net operating range 



- 2 









H ^ A ^ 



+ J 



1 + K - H 



(9.10.11) 



where positive values of A correspond to mode P, negative values to mode W. 

9.11 Slipstream Core Radius 

For any arbitrarily assigned values of k and A which satisfy the constraint 
(9.10.11), parameter Q is fixed by (9. 10. 3) or (9.10.4) and parameter y by 
(9.10.2). Eq. (9.5.6) which fixes the relative radius T|q of the non-swirling 
core in the slipstream can then be non -dimens ionali zed in the form 



80 



(9.1X.1) 




9-12 Parametric Solution for Conventional Performance Coefficients 

With all of the foregoing details fixed, it is a routine matter to non- 
dimens ionalize and integrate the fundamental thrust and power relations (9-6.1) 
and (9-7.1). Upon utilizing relations (9-3) > we finally obtain the following 
parametric solutions for C^, Cp and J. 



As explained above, the two auxiliary parameters y ^.nd which occ\ir 
in these equations are themselves known functions of the two independent 
parameters n and A. Hence in principle these equations suffice to fix a family 
of curves of Cp versus J, with Cp as the parameter which changes from one curve 
to the next. 

9.13 Static Thrust and Power Coefficients 

The reader can easily verify that ■under static conditions (V = O), the 
foregoing parametric solution reduces to the following form. 





(9.12.1) 





(9.12.3) 
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H = J = 0 



0 ^ A ^ 1 



(9.13.1) 

(9.13.2) 



= 1 - </l^ 



(9.13.3) 



— - 1 

W 



(9.13.1^) 



Cp = V ( I ) 



(t? (1 - 1\, j - V 




1 ' 


L ’b J 



+ i T]q^I (9.13.5) 



= 



n 



p - T 



A V (1 - 



2/3 



(9.13.6) 



Values of , Cp and Cp/Cp as calculated from these equations were 
previously presented in Fig. 2.9.I 

9.1U Maximum Attainable Thrust and Power Coefficients 

The reader can also readily verify that under the absolute maximum loading 
limit of Eq. (9.10.8) the general parametric solution reduces to the following 
simplified results. 



^lAX 



= y^i + - H 



*^MAX ” ^ 



= 2 



\) 



(9.1^. 1) 
(9.1^. 2) 

(9. 1^.3) 
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-J/ 



1 + - 




J = HH 




1 + 



K 



+ H 






(9-1J+.6) 



Eq. (9-l*+*3) shows that under these absolute maxinrum loading conditions 
the non-swirling core occupies exactly half the cross-sectional area of the 
slipstreaml Under normal light loadings, the core is of course very much 
smaller than this . 

Ciirves of C„] and Cp| versus J as computed from these relations 
7 MAX 7 MAX 

were previously presented in Fig. 2.9.2. 

9.15 Relative Flow Angle Downstream of Actuator 

•In studies relating to pitch setting of variable pitch actuators, it is 

frequently useful to be able to compute the theoretical relative flow angle 

B, at radius r, on the downstream side of the actuator disc. For the idealized 
d d 

actuator, this is defined by the expression 



tan 




(9.15.1) 



To relate this to corresponding values in the remote slipstream, we note 

that 
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(2r' ^ 



H + A 



H + ^ 



(9.15.2) 






K + ^ 



K 



V 



(9.15.3) 



Upon combining, rearranging and non-dimens ionalizing we finally obtain 




(9.I5.U) 



The corresponding flow angle p in the remote slipstream itself is given 
by the simpler expression 



tan 



H 




(9.15.5) 



Both of these angles vary as fimctions of Tj . To find an angle which is 
respresentative of the pitch setting in some approximate overall sense, it is 
necessary to use an appropriate mean relative radius T| . For example, we may 
choose as a suitable mean value for this purpose the quantity 

/i + 

For an actuator with a sufficiently large number of blades (high blade 
density) the outlet flow angle at radius will approximately equal the 

outlet blade angle itself at this point. For larger blade spacings, the 
relation between flow angle and blade angle is more complex and lies outside 
the scope of this discussion. 
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By combining relations (9. 15-6) with (9.15.4) and imposing the constraint 



it is possible to find a relation between h and A which corresponds approximately 



substituted into the general parametric solution of section 9-12 then roughly 
simulates the performance of an actuator at fixed pitch setting. 

10. Dynamic Stability of the Idealized Lifting Rotor in Vertical Ascent or 



In this section we consider the response of a lifting rotor to small 
perturbations from equilibrium. The treatment is deliberately simplified to 
forcus attention strictly on the basic phenomena. An ideal actuator is assumed, 



vertical climb or descent, except for the imposed perturbations. At the steady 
reference condition, lift L is equal to weight W. But unlike the weight W 
which is constant, lift L is subject to small perturbations 6L about the 
equilibrium value. Lift L and weight W are the only forces considered in this 
analysis, drag and other incidental forces being deliberately ignored. 

Under these circumstances it is advantageous to substitute W for L as 
dimensional reference parameter. Hence we now work in the p, D, W dimensional 
system. The relevant coefficients are defined as follows 



p, = constant 




(9.15.7) 



to constant pitch setting (at high blade density). The corresponding values 



Descent 



with = 0. We analyze a rotorcraft of gross weight W operating in steady 



C 



P 



( 10 . 1 . 1 ) 





H 



( 10 . 1 . 2 ) 
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(10.1.3) 



c 



3L 



L 

W 



In this system, the fundamental equations for the ideal actuator take 
the following form - 
For Modes P, R and W - 




2 



+ a. 
1 




( 10 . 2 . 1 ) 



For Mode T 



C 





- "3L "3H 



( 10 . 2 . 2 ) 



ac. 



For the stability analysis it is necessary to find the partial derivatives 



3P\ /^3P\ 

and 1. This is easily accomplished by differentiating Eqs. (10.2), 

3h/ 



However, in the end we require only the eqviilibrium values of these various 
quantities. The equilibrium values are obtained simply by setting C_ = 1 
in the various expressions. Of course, this may be done only after the 
required differentiations are carried out. In this way the following resxilts 
are obtained 



For modes P, W and R - 



C 



3P 





X = 



3C 



3P 



ac 



3H 







(10.3.1) 



(10.3.2) 
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(10.3.3) 



For mode T - 



C' 



HI +3t 



V _ . . 3P _ H , L\ 2 / TT CO- 
■ dc_ ~ 2 — - 



3L 



H ^ 2 T 

^i\/— / 



n cp 






(10.4.1) 



X - _ JL T 0 ^ 1 



(10.4.2) 



^3P 

Y = — — = - C • 
^3L 



(10.4.3) 



The basic equations (10.2) can be represented symbolically in the form 



C 



3P 





(10.5.1) 



Consequently, the small changes occurring during a perturbation may be written 




or 



6 C 



3P 



= X 6 C 3 ^ H- Y 6 C 3 L 



(10.5.2) 



The equation of motion of the rotorcraft during a small pertiirbation is - 



W 






( 10 . 6 . 1 ) 
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By utilizing the definitions (10.1.2) and (10,1.3) this may be put into the 
form ' 



6C 



3L 



1 _d_ 

g dt 




( 10 . 6 . 2 ) 



It is now advantageous to introduce a non-dimensional time 9 defined as follows 




(10.6.3) 



When expressed in terms of this dimensionless time coordinate 0, the equation 
of motion (10.6.2) becomes simply 



6C 



3L 




(10.6.U) 



It is instructive to include in the analysis the possibility of incorporating 
an automatic power control system. Specifically, we assume that power input 
or output may be modulated in response to pertiirbation velocity and/or accel- 
eration as sensed by appropriate instrumentation. Thus let 



“3P “ - “ d? 


^^3H 


- V 


^3H~ 




_ - 




_ 



(10.7.1) 



where and v are two non-negative constants representing the response of the 
control system to perturbation acceleration and perturbation velocity, 
respectively. Of course, we may set [j, = v = 0, thus simulating a control 
system which maintains constant power independently of the state of the motion 
of the craft. 

Upon substituting Eqs , (10.6.4) and (10.7.1) into (10.5.2), we obtain the 
basic response equation 



88 



d 




/x + v) 
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de 


^^3H 




1 

CO 



The solution of Eq. (9. 8.1) can be taken in the form 



( 10 . 8 . 1 ) 



6C 



3H 




-'0 



( 10 . 8 . 2 ) 



_d_ 

d0 





(10.8.3) 



Upon substituting these last two expressions into the basic response 
equation (10.8.2), we finally obtain the result 



e = - 



X + v) 
Y + ij, 



(10.8.4) 



where the functions X and Y are as defined in Eqs. (10. 3) and (10.4), while jj, 
and V are two non-negative constants which characterize the control system, 
if any. Under conditions of constant power input or output, = v = 0, and 
Eq. (10.8.4) reduces fiirther to 




(10.8.5) 



Stability of operation in vertical ascent or descent requires that p 

be negative. We shall say that Eq. (10.8.5) defines the inherent stability 

of the idealized rotor, while (10.8.4) defines the augmented stability. 

Study of the actuator derivative functions X and Y as given earlier in 

Fig. 2.10.1 shows that X is positive for all positive and negative values of 

the vertical velocity C.,.,. The function Y, on the other hand, turns out to 

3H 
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be positive for values of algebraically greater than - . 1.693, and 

negative for values of C-^v algebraically less than - — . This cross-over 

3“ ^ 

point coincides with point M on the operating envelope. 

We conclude that the ideal lifting rotor is inherently stable over the 
range 

- -Ir ^ C-A s: + 00 (10.8.6) 

This includes all rates of ascent, as well as descent over the full range of 
power-on reversed thrust operation, and descent over an initial small portion 
of the windmilling range. It is to be expected that for any real rotor, the 
range of inherent stability is likely to be considerably more restricted than 
this . 

Eq. (10.8.4) shows very clearly that the range of stable operation can be 
extended as far below the operating point M as desired simply by introducing 
the requisite amount of acceleration response as denoted by the constant |i. 

Note that the velocity response factor v does not affect the range of stability 
but merely increases the rate at which disturbances either grow or decay. 

Typical stability response curves for various combinations of control 
parameters as computed from the foregoing equations were shown earlier in 
Fig. 2.10.2. 

11. Actuator Efficiency and Effectiveness 

In analyzing the performance of a propeller, windmill or lifting rotor, 
we are concerned not only with determining the absolute values of force, power 
and so on, but also with assessing the comparative performance of the device 
in relation to its theoretical potential. In this connection, an essential 
first step is to reduce all performance parameters to appropriate dimensionless 
form since this normalizes the performance of devices of every size and capacity 
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to a common basis of* comparison. However, this first step, while essential, 
is not necessarily sufficient in itself. For example, to know the actual 
values of the conventional coefficients C„, C_ and J in a specific instance 

r r 

is indispensable for some purposes, but nevertheless does not directly indicate 
the degree of departure from ideal performance. Fortunately, this significant 
aspect of actuator performance can be expressed simply and consistently by the 
concept of power loss as defined elsewhere in this paper. Naturally, this 
loss itself should be expressed in dimensionless form by the coefficient C^, 

Cif, or as appropriate, according to the context. 

In some cases it is useful to break the total shaft power P down in some- 

I 

what more detail as shown, for example, in Table 11 . 1 . From this breakdown of 
total shaft power P into the four basic components shown, it is possible to 
define a sequence of significant power ratios as indicated by Eqs. (U.l.l) 
through (11.2.3) of the table. Note that these definitions all employ the 
sign changing operator o^ whose definition is repeated in Eq. (ll.U.l). The 
purpose, of course, is to keep the efficiency and effectiveness ratios as 
defined in the table always less than unity, regardless of the algebraic sign 
of the ideal shaft power P^. 

The key relations summarized in Eqs. (II.3) follow at once from the earlier 
definitions. Notice from Eq. (II.3.3) that the overall power effectiveness ep 
is closely related to the power loss coefficient 

The familiar concept of efficiency, when applied to actuators, is very 
useful in some contexts, but completely useless in others. For example, for 
a propeller operating under cruise conditions, overall efficiency 'll is ordinarily 
a factor of prime importance. On the other hand, for a lifting rotor in 
stationary hovering fli^t, both T| and T|^ vanish identically and obviously fail 
to provide relevant figures of merit. Nevertheless, the power effectiveness 
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factors 6j^, and continue to be significant in this case and do provide’; 
useful figvires of merit. ; 

Notice that the definition of overall power effectiveness ep amounts tol, 
a comparison of ideal and actual power under conditions of equal thrust. Inj^ 
fact the definition (11.2. 3) may be rewritten in equivalent form as 




( 11 . 5 ) 



where the subscript signifies that ideal and actual power coefficients a: 

m I 

compared on the basis of equal force coefficient. ; 

1 

It is equally possible, however, to compare the ideal and actual force ' 
coefficients on the basis of equal power coefficients for both, and in some 
situations this is a preferable method of comparison. We recall that for , 

I 

any assigned value of C^p, the ideal actuator has two distinct values of 

thrust coefficients which we here denote as C.,p\ and C..J\ . At the same 

value of C-_, the actual thrust coefficient C.„ always his within the range 
±r -Lb 






MIN 






IF 



MAX 



(ll.6.]j 



Consequently, there are two significant force ratios which we define as folic, 
namely. 



and 



®F = 



'IF 



'IF 



MAX 



'IP 






'IF 



'IF 



MIN 



-<C 



IP 



( 11 . 6 . 2 ^ 



( 11 . 6.3 
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Notice that the sign changing exponent a is needed in Eq. (11.6.2) to keep e„ 



always less than unity regardless of the sign of C^p, but that no sign changing 



has the same sign as but that is always negative. We may term 

ep the thrust effectiveness for positive values of C^p, or the drag effective- 
ness for negative values of 0^,^. We may term e^. the braking effectiveness in 
all cases. This terminology, while arbitrary, is reasonably convenient and 

clear. The values of e„ or e can be readily inferred in any given case from 

r 

the basic relations presented elsewhere in this paper. 

The power loss coefficients, the efficiencies, or the p>ower and force 
effectiveness ratios defined in this section provide suitable figures of merit 
for virtually every situation. 



exponent is needed in Eq. (11.6. 3). In this connection, .recall that 




MAX 
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Table 11 ♦! Significant Power Components and Ratios 



Components of Total Shaft Power 
Component 
Thrust Power 
Basic Slipstream Loss 
Idealized Rotation Loss 
Residual Loss 



Sign 


Symbol 


Subtotals 
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Total Shaft Power 

Significant Power Ratios 
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P. 
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a. 



o. 

1 



o. 

1 



p. 

1 



Key Relations 



+ P 



Overall Efficiency 



Ideal Efficiency 



Rotation Power Effectiveness 



Residual Power Effectiveness 



Overall Power Effectiveness 
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